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Using data from FOCUS (E831) experiment at Fermilab, we present a model independent partial- 
wave analysis of the K~n + S-wave amplitude from the decay D + — > K~ir + n + . The S-wave is a 
generic complex function to be determined directly from the data fit. The P- and D-waves are 
parameterized by a sum of Breit-Wigner amplitudes. The measurement of the S-wave amplitude 
covers the whole elastic range of the K~n + system. 

PACS numbers: 13.25Ft,13.30Eg,13.87Fh 



I. INTRODUCTION 

Over forty years have passed since the birth of the Constituent Quark Model, yet the scalar mesons still challenge 
theoreticians and experimentalists. Many states have been reported. Some still need confirmation, others need to 
have better measurements of the pole position and couplings to specific channels. The identification of the nature of 
each state — regular qq mesons, tetraquarks, molecules, glueballs — is a major task which will only be accomplished 
combining results from different types of data. 

An important problem is the understanding of the low energy part of the S-wave K~n + spectrum, where the 
existence of an I = 1/2 state, the k(800) meson, has been the subject of a long-standing debate. Evidence for a 
neutral low mass scalar state in heavy flavor decays has been reported by several experiments^], H, EJ B| ■ The pole 
position has been determined recently using Roy-Steiner representations of K~ir + scattering [5]. However, evidence 
for the charged partner is still scarce and conflicting 0, 0] ■ 

The primary source for the Kir — > Kir scattering has been the data from the classic LASS experiment Q, K~p — > 
K~ir + n. With a cut at low momentum transfers, the K~p interaction is assumed to be entirely due to the one-pion- 
exchange mechanism. The incident pion is, therefore, not a real, asymptotically free particle, but a nearly on-shell 
virtual state. An additonal cut on the -K + n mass was set to avoid baryonic intermediate states. The LASS analysis 
was performed on a sample containing 151 thousand events. With this sample LASS found that the Ktt cross section 
is elastic up to the K threshold (1.454 GeV/c 2 ). Unfortunately, LASS data start only at thk-k = 825 MeV/c 2 . 

Heavy flavor decays are currently the only way to access the whole elastic range of the Kir spectrum, starting from 
threshold. A golden mode for the neutral Kir system is the decay D + — > A'~7r + 7r + , which has a largely dominant 
S-wave component — a common feature of three-body final states with identical pions. 

This decay was already studied by the Fermilab FOCUS Collaboration Q. In our previous work the D + — > K~ii + it + 
Dalitz plot was analyzed with the K-matrix formalism, which was applied for the first time in Dalitz plot analysis of 
D decays in the FOCUS study of the D + — ► tt + it~it + decay As a cross check, a fit with the usual isobar model 
was also performed. 

In the isobar model the S-wave is represented by a coherent sum of a uniform nonresonant term plus two relativistic 
Breit-Wigner amplitudes. A good fit can be achieved, but it is difficult to determine the relative amount of each 
S-wave component. In order to illustrate the correlation between the S-wave components of the isobar model an 
ensemble of 2,000 D + — >■ K~ir + ir + Dalitz plots was simulated using the set of parameters from our isobar fit (Table 
2 of reference [3(). Each simulated Dalitz plot was fitted and a scatter plot of the nonresonant versus k(800)7t + decay 
fractions is presented in Fig. [TJ One can clearly see that a better description of the S-wave requires one to go beyond 
the isobar model. 

The K-matrix formalism is based on the assumption that there is no three-body final state interaction. In this 
approach the dynamics of the K~tt + it + final state are driven by the K~ir + system. Data on D decays and on 
scattering would be directly related and the two-body unitarity would become a constraint. The evolution of the 
K~ir + pair is fixed to that of Ktt clastic scattering, considering the contribution of both 1=1/2 and 7=3/2 K~ir + 
amplitudes. The parameters of the K~iv + production amplitude and the relative amount and phase of the two isospin 
components are determined by the fit. A good description of the data was obtained, with an important contribution 
of the 7 = 3/2 Ktt amplitude. The production amplitude has a slowly varying phase. The conclusion of this study is 
that data on D decays and scattering are consistent. The three-body final state interactions would, therefore, play a 
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FIG. 1: Correlation between nonresonant and nn + decay fractions in the isobar model. An ensemble of D + — > K~n + n + 
Dalitz plots was simulated and fitted with the isobar model. The scatter plot shows a high correlation between the two largest 
contributions of the S-wave. 



marginal role. 

In this paper we complete our study of the D + — > K~ir + ir + (charge conjugate states are always implied) Dalitz 
plot, applying, to the same data set, the model-independent partial wave analysis (MIPWA) technique, developed by 
the E791 Collaboration 0. In this method, the K 7r + S-wave amplitude is parameterized by a generic complex 
function, to be determined directly from the data. The only assumption common to all other Dalitz plot analyses is 
that the P- and D-waves are well represented by a sum of Breit-Wigner amplitudes. The K~tt + spectrum is divided 
into slices. The magnitude and phase of the S-wave component at the edge of each slice are determined by the fit. A 
cubic spline interpolation is used to obtain the S-wave magnitude and phase at any point in the spectrum. 

The MIPWA technique provides a model-independent way to determine the K~tt + S-wave amplitude. The result, 
however, is inclusive. The measured phase, in addition to the 1=1/2 K~ir + phase, may contain contributions from 
the 7=3/2 components, as well as possible contributions from three-body final state interactions. 

The paper is organized as follows. In the next section we describe the selection of the data sample. The MIPWA 
formalism is described in Section III. The results of the MIPWA fit are presented in Section IV. 



II. THE D + -> K~n + n + SAMPLE 



FOCUS is a charm photo-production experiment which collected data during the 1996-97 fixed target run at 
Fermilab. The photon beam was produced by means of bremsstrahlung, from electron and positron beams (typically 
with 300 GeV endpoint energy). The electron/positron beams were obtained from the 800 GeV Tevatron proton 
beam. The photon beam interacted with a segmented BeO target [ill ]. The mean photon energy for reconstructed 
charm events is ~ 180 GeV. 

The FOCUS spectrometer has a system of three multi-cell threshold Cerenkov counters to perform the charged 
particle identification, separating kaons from pions up to a momentum of 60 GeV/c. The identification and separation 
of charm primary (production) and secondary (decay) vertecis are made by two systems of silicon micro-vertex 
detectors. The first system consists of 4 planes of micro-strips interleaved with the experimental target [11] and the 
second system consists of 12 planes of micro-strips located downstream of the target. The charged particle momentum 
is determined by measuring the deflections in two magnets of opposite polarity through five stations of multi-wire 
proportional chambers. 

The data set used in this analysis is the same as in Ref. [3|. The final states are selected using a candidate driven 
vertex algorithm [l3l ] . A secondary vertex is formed from the three candidate tracks. The momentum of the resultant 
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FIG. 2: The K tt + tv + invariant mass distribution. 



D + candidate is used as a seed track to intersect the other reconstructed tracks and to search for a primary vertex. 
The primary vertex must have least two reconstructed tracks in addition to the D + seed. The confidence level of 
each vertex is required to be greater than 1%. Once the production and decay vertecis are determined, the distance L 
between the vertecis and its error ox are computed. The quantity L / gl is an unbiased measure of the significance of 
detachment between the primary and secondary vertecis. This is the most important criterium for separating charm 
events from non-charm prompt backgrounds. Signal quality is further enhanced by isolation requirements. Tracks 
forming the D candidate vertex must have a confidence level smaller than 0.001% to form a vertex with the tracks 
from the primary vertex. In addition, all remaining tracks not assigned to either the primary or the secondary vertex 
must have a confidence level smaller than 0.1% to form a vertex with the D candidate daughters. 

Particle identification cuts used in FOCUS are based on likelihood ratios between the various particle identification 
hypotheses. These likelihoods are computed for a given track from the observed firing response (on or off) of all the 
cells that are within the track's (j3 = 1) Cerenkov cone for each of our three Cerenkov counters. The product of all firing 
probabilities for all the cells within the three Cerenkov cones produces a x 2 -like variable Wi — —2 ln(Likelihood) where 
% ranges over the electron, pion, kaon and proton hypotheses • The kaon track is required to have Ak — W^ — Wk 
greater than 3; both pion candidates are required to satisfy = Wk — W^ greater than 3; in addition, all tracks 
are required to be separated by less than 5 units from the best hypothesis, that is AW = W m in — Wk < 5 and 
AW = W m in — W-n < 5 . These Cerenkov cuts reduce the contamination of Df — > K~K + ir + background to a 
negligible level. 

Using the set of selection cuts just described, we obtain the invariant K~7t + n + mass distribution shown in Fig. 
[5J The mass plot of Fig. [5] is fitted with a function that includes two Gaussian functions with different widths 
and the same mean, which take into account differences in the resolution in the momentum determination of our 
spectrometer and an exponential function for the background. The events used in the MIPWA fit correspond 
to the shaded area in Fig. [21 i.e., events with 1.8515 < Mkttit < 1.9031 GeV/c 2 . Events in this mass region that 
lie outside the kinematic limit defined by the nominal D + mass are discarded. The final data subset contains 53,595 
events, with a purity (S/(S+B)) of 98.8%. 

The symmetrized Dalitz plot of these events (two entries per event) is shown in Fig. [3j A narrow band corresponding 
to the D + — > A"*(892)°7T + events can be clearly seen. The asymmetry in each K* (892)° lobe is evident and it is 
caused by the interference between this state and the K~tt + S-wave. Indeed, it is this interference with the P-wave 
that allows one to access the S-wave phase. 
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FIG. 3: The K~n+n+ Dalitz plot. 



III. THE MODEL INDEPENDENT PARTIAL WAVE ANALYSIS FORMALISM 

In the MIPWA formalism the Dalitz plot of Fig. [3]is described by a coherent sum of three partial waves, correspond- 
ing to the K~ir + system in the angular momentum states L =0, 1 and 2. The partial waves are complex functions 
of the two Kir invariant masses squared, s a = {px + Pir a ) 2 and s b = (pk + P-n b ) 2 i which specify the kinematics of the 
D + — ► K^ir+ir^ decay. Each partial wave is Bose-symmetrized with respect to the identical pions, 



A L = A L (s a ,s b ) + A L (s b , s a ). 
The K ~tt + S-wave amplitude is an unknown complex function of the K ~7r + mass squared, 



(1) 



A (s a , s b ) = a Q ( Sa )e^ is ^ + a (s b )e l ^l 



(2) 



No assumption about the content of the S-wave is made: the real functions ao(s) and (f>o(s) are determined directly 
by the Dalitz plot fit. The K~tt + mass spectrum is divided into 39 slices of the same size. For each of the 40 
endpoints Sk there are two free parameters, a* and <pk, defining the function Ao(s a ,s b ) at that position. A cubic 
spline interpolation is used to define the values of both ao(s) and 4>o(s) between Sk < s < Sk+i- The S-wave has, 
therefore, a set of 40 pairs (a^, </>&) of fit parameters. 

The K~tt + P-wave amplitude has two components, namely the K (892)°7r + , taken as the reference mode, and the 
iT(1680) 7r+, 



Ai{s a , Si) = Ff (s a , s b )F^(s a , s b )[c BW K *(892)(sa) + ciBW K *( lfim) (s a )]Mi{s a , s b ) 
The D-wave has only one component, the i^ 2 (1430)7r + mode, 



(3) 



A 2 (s a ,s b ) =c 2 [F°(s a , S b )F?(Sa, s b )BW K > {li30) (s a )] M 2 (s a: s b ) 



(4) 



The complex coefficients Ci are also fit parameters, except for Cq, the coefficient of the reference mode, which is 
fixed to 1.0. 

In the above equations F D and F R are the usual Blatt-Weisskopf form factors [l5| , 
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F L=1 = [l + (rqf]- 1/2 , ( 5 ) 

and 

F L = 2 = [9 + 3(rq) 2 + (rq) 4 )]- 1 / 2 . ( 6 ) 

where q is the momentum of the resonance decay products in the resonance rest frame. The form factor parameters 
r = Td for the D decay vertex, and r = Tr for the resonance decay, are fixed at the values used in Ref. 
r D = 1.5 (GeV/c)" 1 and r R = 5.0 (GeV/c)" 1 . 

The functions M. l are the spin amplitudes, accounting for angular momentum conservation. For a spin-1 resonance 
D + — ► Rn^ , R — > K~tt+ the corresponding spin amplitude is 

M L =i = Y,P%MP> M ) e «(P,M)(p na -p K y, (7) 

M 

where e M is the resonance polarization vector with magnetic quantum number M, p is the momentum 4- vector and 
pf are the momenta of the final state particles. After summing over the unobserved resonance polarization states, 
the spin amplitude reduces to 

Ml=i = -2 | p-K h || Pk | cos 9, (8) 

where 9 is the cosine of the angle formed by px and p 7!h in the resonance rest frame. 
In the case of the if 2 (1430)7r + mode, the spin amplitude is 

M L =2 = Y J P^^ b e ^ M ) e ^ M ) e ^M)e {p>,M)(p^ a -p K ) a (pn a ~p K f, (9) 

M 

which, after summing over the resonance polarization, reduces to 

M L=2 - |(| p W6 ||pjr |) 2 (3cos 2 - 1). (10) 
The relativistic Breit-Wigner has an energy dependent width, 

BW= -rVr (V ( U ) 

S - S Q + ly/So 1 tot(s) 

where s is the K~it + mass squared, sq the resonance nominal mass and 

'so f 1 \ F l(<l) 



s 



r tot ( s ) = r W- M- (12) 



where L is the orbital angular momentum in the rest frame of the decaying resonance. 

The signal distribution is corrected on an event-by-event basis for the acceptance. The acceptance includes geometry, 
detector and selection cuts efficiency. It is determined by a full Monte Carlo simulation of events: the j—N interaction, 
event propagation through the spectrometer, event reconstruction and selection of the Kitir sample. The acceptance 
function is obtained by fitting the Dalitz plot of Monte Carlo events to a 10th order polynomial. 

The signal probability distribution is normalized to unity, 

1 I .2 

Ps(s a ,s b ) = —e(s a) s b ) Y]Al , (13) 
i\ s I 

where s(s a , Sb) is the acceptance function and Ns the overall normalization constant, 



7 



N s = I ds a ds b e{s a ,s b )\y^A L \ ■ (14) 



The background probability distribution is fixed in the fit. The background shape is determined by a fit to the 
Dalitz plot of events from the K + tt~tt~ mass sidebands Q. The signal fraction is estimated by a fit to the K + ir~Tr~ 
mass spectrum. 

In the MIPWA fit there are 40x2 + 2x2=84 free parameters. The optimum set of parameters is determined by an 
unbinned maximum likelihood fit, minimizing the quantity w = — 21n(L), where the likelihood function, L, is given 
by 

L= J] [fsP l s(Sa,S b ) + (l-f S )Ph(Sa,S b )] (15) 
events 

where f s is the signal fraction f s = S/(S + B). 

Decay fractions are obtained from the coefficients Ck, determined by the fit, and after integrating the overall signal 
amplitude over the phase space, 

, J ds a ds b \c k Ak(s a , s b )\ 2 

Jk — 2 ■ V ib / 



/ ds a ds b Y^j c jAj(s a ,s b ) 



Errors on the fractions include errors on both magnitudes and phases, and are computed using the full covariance 
matrix. 



IV. RESULTS OF THE MIPWA 



The decay fractions resulting from the MIPWA fit are presented in Table [IJ For comparison, the third and fourth 
columns have the fractions from our previous fits using the K-matrix formalism and the isobar model (l7l |. 

The data is well described by a P-wave model with two components. No improvement in the fit quality is observed 
when a third component, the K (1410)7r + mode, is added. The contribution of this mode is consistent with zero. 
As in previous analyses of the D + — > K + t:~tt~ Dalitz plot, the S-wave component is dominant. The decay fractions 
from the MIPWA and from our previous K-matrix Dalitz plot fit are in good agreement. 

In Table HT1 the MIPWA decay fractions are compared to the ones from E791 and CLEO-c. Our results are in good 
agreement with the decay fractions from E791. The total S-wave contribution from CLEO-c is significantly higher if 
we add the binned and the X (1430)7r + fractions. 

The fitted values of the S-wave magnitudes and phases are presented in Table [TTT1 and plotted in Figs. 2] and 03 The 
error bars in these figures contain the statistical and systematic errors folded in quadrature. The dashed line in Fig. 
S] indicates the Krf threshold, the upper limit of the region where the K~tt + amplitude is predominantly elastic 

The S-wave phase grows continuously across the elastic region, starting at -138° and with a total variation of 200°. 
After a sudden drop near the K (U30) mass, the phase becomes nearly constant. 

The S-wave magnitude is a decreasing function up to mjfj ~ 1.2 GeV/c 2 . There is a dip near the K (U30) mass, 
which is most readily explained by the interference between the different components of the S-wave. 

The measured magnitudes are more affected by the systematic uncertainties than are the measured phases. In both 
cases the systematic uncertainties are comparable to or larger than the statistical errors. 



A. Goodness-of-fit 



For all fits the goodness-of-fit is accessed through a two-dimensional \ 2 test, using an adaptive binning algorithm. 
The folded Dalitz plot is divided into 844 cells of variable size, with a minimum occupancy of 50 data events, in such 
a way that all cells have a nearly equal and sufficiently large population. This procedure allows us to test the fit 
quality in great detail across the Dalitz plot. For each cell we define the x 2 as 



2 (^-obs ^exp) 



exp 



(17) 
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TABLE I: Decay fractions (%) from the MIPWA Dalitz plot fit. In the MIPWA column the first error is statistical, the second 
and third errors are, respectively, our estimate of the split sample and fit variant systematic uncertainties, and the last error is 
the systematic error due to the uncertainty in the parameters of the other waves. Below the decay fractions and the phases, in 
degrees, are indicated. 



mode 



FOCUS MIPWA 



FOCUS K-matrix 



FOCUS isobar model 



K 7r S-wave 
Z*(892) 7r+ 

¥'(1410)%+ 

7T(1680) 7r + 

^(1430)7t+ 



80.24±1.38±0.23±0.25±0.26 

12.36±0.34±0.19±0.16±0.23 
(fixed) 



1.75±0.62±0.24±0.23±0.42 
67±6±2±2±3 

0.58±0.1±0.04±0.03±0.04 

336±7±3±2±2 



83.23±1.50±0.04±0.07 

13.61±0.41±0.01±0.30 
(fixed) 

0.48±0.21±0.012±0.17 
293±17±0.4±7 

1.90±0.65±0.01±0.43 
1±7±0.2±6 

0.39±0.1±0.004±0.05 
296±7±0.3±1 



13.7±0.4±0.6±0.3 
(fixed) 

0.2±0.1±0.I±0.04 
350±34±17±15 

1.8±0.4±0.2±0.3 

3±7±4±8 

0.4±0.05±0.04±0.03 
3I9±8±2±2 

17.5±1.5±0.8±0.4 
36±5±2±1.2 



22.4±3.7±1.2±1.5 
199±6±1±5 



nonresonant 



29.7±4.5±1.5±2.1 
325±4±2±1.2 



TABLE II: Decay fractions (%) and phases, in degress, from the MIPWA Dalitz plot fit compared to E791 and CLEO- 



mode 



FOCUS MIPWA 



E791 



CLEO-c 



"/? 7r + S-wave 
Zo(1430)tt + 

Z*(892)°7r+ 

7T(1680)°7r + 

7?;(1430)7r+ 



.24±1.38±0.23±0.25±0.26 



12.36±0.34±0.19±0.16±0.23 
(fixed) 

I.75±0.62±0.24±0.23±0.42 

67±6±2±2±3 

0.58±0.1±0.04±0.03±0.04 

336±7±3±2±2 



78.6±2.3 



11.9±2.0 

(fixed) 

1.2±1.2 
43±17 

0.2±0.1 
-I2±29 



83.8±3.8 

13.3±0.62 
51 (fixed) 

9.88±0.46 
(fixed) 

0.20±0.12 
113±14 

0.20±0.04 

15±9 



In the above expression n eX p is the expected population of each cell, given by a Monte Carlo simulation performed 
with 1,000,000 events generated according to the model resulting from the MIPWA fit, and cr oxp is the uncertainty 
on this number. The overall x 2 is a sum of the xf over all cells. The number of degrees-of-freedom is given by the 
number of cells minus the number of fit parameters. From these two quantities we estimate the confidence level of 
our fits. 

The overall x 2 of the MIPWA fit is x 2 =818.8 (844-84=760 degrees of freedom), which corresponds to a confidence 
level of 6.8%. The x 2 distribution across the Dalitz plot is shown in Fig. [6] 

The Dalitz plot projections (highest and lowest K~tt + invariant mass squared) are plotted in Fig. [7j with the fit 
result superimposed (solid histograms). 
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TABLE III: Magnitudes and phases of the S-wave from MIPWA fit. The first error is statistical. The second and third errors 
are, respectively, our estimate of the split sample and fit variant systematic uncertainties, whereas the last error is the systematic 
error due to the uncertanty in the parameters of the other waves. The full systematic error is a sum in quadrature of these 
three errors. This is the number between parentheses. 
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B. Systematic uncertainties 

Systematic uncertainties may come from different sources. We have performed split sample studies, in which the 
data was divided into four sets of independent samples, according to the parent D meson charge and momentum. The 
split sample component takes into account the possible systematics introduced by a residual difference between data 
and Monte Carlo, due to a possible mismatch in the reproduction of the D + production. A technique, employed in 
FOCUS and modeled after the S-f actor method from the Particle Data Group [l6j], was used to try to separate true 
systematic variations from statistical fluctuations. We found a small effect from the split sample studies. 

A second class of studies is the fit variant, in which the fit of the whole data set is performed under different 
conditions. Fit variants included changes in the background level and in the first derivatives of the spline at the edges 
of the K~7r + spectrum. The fit variant component can be estimated by the r.m.s. of the measurements. 

The third and dominant source of systematic errors comes from the uncertainty in the parameters of the P- and 
D-waves. This includes uncertainties on the values of the parameters tr and Td- We repeated the fit changing by 
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FIG. 4: The S-wave phase as a function of the K~n + mass from the MIPWA K~~ir + ir + Dalitz plot fit. The hashed vertical 
line shows the elastic range according to LASS. 
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FIG. 5: The S-wave magnitude as a function of the K -k + mass from the MIPWA K -k + tv + Dalitz plot fit. 



±1(7, one at a time, the values of the mass and width of the high-mass vector resonances, according to the PDG, and 
of the parameters rn and rrj>. This component is also estimated by the r.m.s. of the measurements. 

The contributions of each source are quoted individually in Table HI |TT] and IIIIl The overall systematic uncertainty 
was obtained adding in quadrature the three components described above, and corresponds to the values in parentheses 
in Tabic [EE 
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FIG. 6: The x 2 distribution across the folded Dalitz plot. 




FIG. 7: Comparison between the K tt + tt + Dalitz plot projections and the MIPWA fit, with the lowest (left plot) and highest 
(right plot) K~ir + invariant mass squared. In the plots the solid histogram is a projection of the fit. 

V. SUMMARY AND CONCLUSIONS 

A Dalitz plot fit was performed with the MIPWA technique. The K~ir + S-wave amplitude was determined directly 
from data, with no assumption about its nature. The only hypotheses are that the decay amplitude can be described 
by a sum of partial waves, and that the P- and D- waves are well described by a coherent sum of Breit-Wigner 
amplitudes. 

The MIPWA decay fractions are in good agreement with our previous analysis and with the E791 results. A large 
dominance of the S-wave component is observed in this decay. 

The phase of the S-wave amplitude grows continuously across the elastic range, with a total variation of approxi- 
mately 200°. At the K~tt + threshold there is a phase difference of approximately -140° between the S- and P-waves. 

The phase variation of the S-wave measured in this analysis and that of E791 agree well, specially in the elastic 
range. Our definition of the S-wave amplitude, eq. [2l differs from that of E791. The latter includes a Gaussian form 
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FIG. 8: Comparison between the S-wave phase from the three different FOCUS fits of the K~ir ir Dalitz plot. Points with 
error bars are the result of the MIPWA fit. The solid line is the central value of the isobar fit. The dashed line is the result of 
the K-matrix fit. 



factor, so one should compare the S-wave magnitude from our analysis to the product of the E791 Gaussian form 
factors and magnitude. We also find a qualitative agreement between the S-wave magnitude measured by the two 
experiments. 

FOCUS has performed a comprehensive study of the D + — > K~ix + n + Dalitz plot. Using the same events, fits with 
the isobar model, the K-matrix formalism and the MIPWA were performed. The three fits have equivalent goodness- 
of-fit. The decay fractions from all fits are in good agreement. In the isobar model there is a strong correlation 
between the nonresonant and kit modes. Although a good fit with this model is achieved, it is difficult to disentangle 
the contribution of these two modes. 

In Fig. [8] the S-wave phase from the three fits are compared. All fits show a good agreement in the interval 
1 < < 1-35 GeV/c 2 . The MIPWA phase is lower than those from the isobar/K-matrix fits for rriK-n < 1 GeV/c 2 . 

In the high mass region the rapid variation of the phase is more pronounced in the isobar/K-matrix fits than in the 
MIPWA. 

The S-wave magnitude from the three fits are compared in Fig. [5] In the isobar and K-matrix fits there is a 
broad maximum at around 0.9 GeV/c 2 , which is absent in the MIPWA fit. In the region 1.2 < rajc* < 1-4 GeV/c 2 
the MIPWA magnitude has a bump whereas in the isobar and K-matrix the magnitude decreases. In the high mass 
region, after the minimum, the magnitude from the MIPWA fit has a steeper variation than that of the isobar and 
K-matrix. 

The D + — > K~ir + ir + decay offers an opportunity to access the K~n + S-wave amplitude near threshold. Except for 
heavy flavor decays, no new data on the K~ir + system are foreseen. The ultimate goal is to extract the 1—1/2 K~ir + 
elastic amplitude, where all resonances are contained. The result of the MIPWA fit, however, may include other effects, 
such as a possible contribution of the 7=3/2 amplitude, or an energy dependent phase introduced by three-body final 
state interactions. The road from the MIPWA S-wave to the 1—1/2 K~ir + elastic amplitude is, unfortunately, not 
direct. Input from theory is necessary. At this level of statistics we are already limited by systematics, which are 
dominated by the uncertainties on resonance parameters. 
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FIG. 9: Comparison between the S-wave magnitude from the three different FOCUS fits of the K~ir tv Dalitz plot. Points 
with error bars are the result of the MIPWA fit. The solid line is the central value of the isobar fit. The dashed line is the 
result of the K-matrix fit. 



the physics departments of the collaborating institutions. This research was supported in part by the U. S. National 
Science Foundation, the U. S. Department of Energy, the Italian Istituto Nazionale di Fisica Nucleare and Ministero 
della Istruzione Universita e Ricerca, the Brazilian Conselho Nacional de Desenvolvimento Cientffico e Tecnologico 
and FAPERJ, CONACyT-Mexico, and the Korea Research Foundation of the Korean Ministry of Education. 



[1] E.M. Aitala et al. (E791 Collaboration), Phys. Rev. Lett. 89, 121801 (2002). 

[2] M. Ablikim et al. (BES Collaboration), Phys. Lett. B598, 149 (2004). 

[3] J.M. Link et al. (FOCUS Collaboration), Phys. Lett. B653, 1 (2007). 

[4] G. Bonvicini et al. (CLEO-c Collaboration), Phys. Rev. D78, 052001 (2008). 

[5] S. Descotes-Genon and B. Moussallam, Eur. Phys. J. C48, 553 (2006). 

[6] B. Aubert et al. (BaBar Collaboration), Phys. Rev. D76, 011102 (2007). 

[7] D. Epifanov et al. (Belle Collaboration) . larXiv:0706.2231 1 

[8] D. Aston et al. (LASS Collaboration), Nucl. Phys. B296, 493 (1988). 

[9] J.M. Link et al. (FOCUS Collaboration), Phys. Lett. B585, 200 (2004). 
[10] E.M. Aitala et al, (E791 Collaboration), Phys. Rev. D73, 032004 (2006). 
[11] P.L. Frabetti et al, (E687 Collaboration), Nucl. Instrum. Meth. A329 (1993) 62. 
[12] J.M. Link et al, (FOCUS Collaboration), Nucl. Instrum. Meth. A516(2004) 364. 
[13] P.L. Frabetti et al, (E687 Collaboration), Nucl. Instrum. Meth. A320 (1992) 519. 
[14] J.M. Link et al, (FOCUS Collaboration), Nucl. Instrum. Meth. A484 (2002) 270. 
[15] J.M. Blatt and V.F. Weisskopf, Theoretical Nuclear Physics (John Wiley & Sons, New York, 1952). 
[16] C. Amsler et al. (Particle Data Group), Phys. Lett. B667, 1 (2008). 

[17] The two statistical errors here reported on the K (892)°7r + decay fraction for both the K-matrix and isobar fits are 
smaller than those quoted in Ref. [j§|. They were overestimated in Ref. [3] as a consequence of a minor mistake in the 
error propagation code, which affected their computation only and not those of the other fit fractions. The two new values 
reported here are the correct ones, and should replace the old ones. 



